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We report on a high statistics measurement of the total and differential cross sections of the 
process 77 — > ■k'^'k~ in the tt^tt^ invariant mass range 0.8 GeV/c? < W < 1.5 GeV/c? with 
85.9 fb~^ of data collected at y/s = 10.58 GeV and 10.52 GeV with the Belle detector. A clear signal 
of the /o (980) resonance is observed in addition to the /2 (1270) resonance. An improved 90% 
confidence level upper limit Z3(?7'(958) ^ tt+tt") < 2.9 x 10"^ is obtained for P- and CP-violating 
decay of the »7'(958) meson using the most conservative assumption about the interference with the 
background. 
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I. INTRODUCTION 



The nature of low mass^mesons remains poorly understood in spite of decades of theoretical 
and experimental effort [jj. In particular, low mass scalar mesons (below 1 GeV/c^) are not 
yet well established experimentally except for the /o (980) and ao (980) mesons, while the 
extensively discussed a (/o(600)) and k (-ft'*(800)) mesons still remain controversial states Q. 
A B factory is well suited for detailed investigations of low mass mesons through two- 
photon production, where overwhelming statistics can be obtained. Two-photon production 
of mesons has advantages over meson production in hadronic processes; the production rate 
can be reliably calculated from QED with r^.y as the only unknown parameter. In addition, a 
meson can be produced alone without additional hadronic debris, and the quantum numbers 
of the final state are restricted to states of charge conjugation C — +1 with J — 1 forbidden 
(Landau- Yang's theorem 

In the past, extensive studies of low mass mesons through 77 — > tttt scattering have 
been made at e+e" colliders: Crystal Ball [i], Mark II 0, JADE Q, TOPAZ MD-1 [8], 
CELLO (ofl and VENUS yLfl]; see Ref. Q for a list of the earlier experiments. Using data from 
Mark II, Crystal Ball, and CELLO, Boglione and Pennington (BP) performed an amplitude 
analysis of 77 tt+tt^ and 77 7r°7r° cross sections fll|. They found two distinct classes 
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of solutions where one solution has a peak ("peak" solution) and the other has a wiggle 
("dip" solution) in the /o (980) mass region. The two solutions give quite different results 
for the two-photon width of the fo (980) and the size of the S-wave component. Thus, it is 
necessary to distinguish them experimentally. 

In this paper, we report on a measurement of the cross sections for the reaction 77 
7r+7r~ with high statistics that are more than two orders of magnitude larger than that of the 
past experiments. The analysis is based on data taken with the Belle detector at the KEKB 
asymmetric-energy (3.5 GeV on 8 GeV) e+e~ colhder The data sample corresponds 

to a total integrated luminosity of 85.9 fb~^, accumulated on the T(4S') resonance {^/s ~ 
10.58 GeV) and 60 MeV below the resonance (8.6 fb~^ of the total). Since the cross section 
difference between the two energies is only about 0.3%, we combine both samples We 
observe the two-photon process e+e^ e+e~7r+7r~ in the "zero-tag" mode, where neither 
the final-state electron nor positron is detected, and the tt+tt" system has small transverse 
momentum. We restrict the virtuality of the incident photons to be small by imposing a 
strict requirement on the transverse-momentum balance of the final-state hadronic system 
with respect to the beam axis. Some of the results reported here are the subject of a separate 
paper focusing on the properties of the /o(980) meson [l^. 

This paper is organized as fohows. A brief description of the detector is given in section ITTl 
The selection criteria are listed in section IIIII There is a well known difficulty in discrim- 
inating /i^ from TT^ in the low momentum region (< 0.8 GeV/c); Section [iVl presents the 
method of particle identification, in particular the method of n/ir separation that we use. 
Evaluation of the detection and trigger efficiencies is described in section |Vl The total and 
differential cross sections are given while their systematic errors are estimated in section IVll 
In section IVIH the resulting spectrum is fitted to obtain the resonance parameters of the 
/o(980) meson and to check consistency in the /2(1270) region. Section rVIIII summarizes 
the results. Appendix [A] gives a detailed description of the background subtraction. Values 
of the total cross sections are given in Appendix [b1 

II. THE BELLE DETECTOR 

The Belle detector is a large-sohd-angle magnetic spectrometer having good momentum 
resolution and particle identification capability in the energy region of interest Here 
we briefiy describe the Belle detector components. Charged track coordinates near the col- 
lision point are measured by a 3-layer sihcon vertex detector (SVD) that surrounds a 2 cm 
radius beryllium beam pipe. Track trajectories are reconstructed in a 50-layer central drift 
chamber (CDC), and momentum measurements are made together with the SVD. An array 
of 1188 sihca-aerogel Cherenkov counters (ACC) provides separation between kaons and 
pions for momenta above 1.2 GeV/c. The time-of-fiight counter (TOF) system consists of 
a barrel of 128 plastic scintillation counters and is effective for K/tt separation for tracks 
with momenta below 1.2 GeV/c. Low energy kaons and protons are also identified through 
specific ionization (dE/dx) measurements in the CDC. Photon detection and energy mea- 
surements of photons and electrons are provided by an electromagnetic calorimeter (ECL) . 
It is comprised of an array of 8736 CsI(Tl) crystals all pointed toward the interaction point, 
greatly enhances the electron identification capability provided through a comparison of 
energy measured in the ECL and momentum in the CDC. These detector components are 
located within a superconducting solenoid coil that provides a uniform magnetic field of 
1.5 T. An iron fiux- return located outside the solenoid coil is instrumented to detect 
mesons and to identify muons (KLM). The z axis of the detector is defined to be opposite to 
the direction of the positron beam. These detector components cover a polar angular range 
between 17° and 150°. 
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III. EVENT SELECTION 



Signal candidates are primarily triggered by a two-track trigger that requires two CDC 
tracks with associated TOF hits and ECL clusters with an opening angle greater than 
135°. Exclusive e+e^ — > e+e^7r+7r^ events are selected by requiring two oppositely charged 
tracks coming from the interaction region; each track is required to satisfy dr < 0.1 cm and 
\dz\ < 2 cm, where dr (dz) is r (z) component of the closest approach to the nominal collision 
point. Here, r is the transverse distance from the z axis. The difference of the dz's of the two 
tracks must satisfy the requirement \dz+ — dz-\ < 1 cm. The event must contain one and 
only one positively charged track that satisfies pt > 0.3 GeV/c and —0.47 < cos9 < 0.82, 
where pt and 9 are the transverse component of momentum and the angle with respect 
to the z-axis. The scalar sum of the track momenta in each event is required to be less 
than 6 GeV/c, and the sum of the ECL energies of the event must be less than 6 GeV. 
Events should not include an extra track with pt > 0.1 GeV/c. The cosine of the opening 
angle of the tracks must be greater than —0.997 to reject cosmic-ray events. The sum of 
transverse momentum vectors of the two tracks {J^Pt) should satisfy \ J2Pt l< 0.1 GeV/c; 
this requirement separates exclusive two-track events from quasi-real two-photon collisions. 



IV. PARTICLE IDENTIFICATION AND fi/n SEPARATION 

Electrons and positrons are clearly distinguished from hadrons using the ratio E/p, where 
E is the energy measured in the ECL, and p is the momentum from the CDC. Kaon (proton) 
candidates are identified using normalized kaon (proton) and pion likeHhood functions {Lk 
(Lp) and L,r, respectively) obtained from the particle identification system (combining the 
information of the CDC, TOF, ACC and ECL) with the criterion Lk/{Lk + L^) > 0.25 
{Lp/{Lp + Ltt) > 0.5), which gives a typical identification efficiency of 90% with a pion 
misidentification probability of 3%. All charged tracks that are not identified as electrons, 
kaons or protons are treated as pions. We require both tracks to be pions. The resulting 
invariant mass (W) distribution is shown in Fig.[TJ The W bin size is chosen to be 5 MeV/c^, 
while the mass resolution is about 2 MeV/c^ according to GEANT-3 [lB| based Monte Carlo 
(MC) simulation. A clear signal corresponding to the /o(980) meson is seen along with the 
well known /2(1270) resonance. 




W (GeV/c') 



FIG. 1: The invariant mass distribution of unseparated (a^ M 4- tt tt^) events. A clear signal for 
the /o (980) can be seen. The large peak around W — 1.2 GeV/c^ corresponds to the /2 (1270) 
resonance. Note the suppressed zero on the vertical scale. 
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In this measurement, the KLM detector cannot be used for muon identification, since it 
is insensitive in the region of interest where the transverse momenta of tracks are below 
0.8 GeV/c. Therefore, we have developed a method for separating tt+tt" and events 
statistically using ECL information; muons deposit energy corresponding to the ionization 
loss for minimum ionizing particles, while pions give a wider energy distribution since they 
may interact hadronically in the ECL, which corresponds to approximately one interaction 
length of material. Typical two-dimensional distributions vs. E_) of the energy deposit 



E± in the ECL for and tt+tt pairs produced by MC are shown in Figs. 2(a) and 2(b) 




FIG. 2: Typical distributions of the energy deposit {E+ vs. E-) in the ECL. Events are produced 
using MC simulation in a range 1.0 GeV/c^ <W < 1.02 GeV/c^ and | cos6i*| < 0.1. 



Probability density functions (PDFs) for the distributions of energy deposits from tt+tt" 
lyjX^ pairs P^^}^- iE+,E-) (-P^+^- (-^+) -£'-)) ^-i^^ obtained with MC simulation. Here i 
represents the i-th bin of {W, |cos0*|) in 20 MeV/c^ and 0.1 steps, where W is the invariant 
mass of the tt+tt^ (or ^^fi") pair in each event (the pion mass is assumed in the calculation), 
and 6* is the polar angle of the produced tt* meson (or lepton) in the center-of-mass 
system of two initial photons. Note that using this method the effect of muons from pion 
decays is taken into account by the pion PDFs. We obtain r^*\ the fraction of in the 

i-th bin through the equation: 

<L (£;+,£_) = iV« (£;+,£-) + (1 - rW)Pit- E^)) , (1) 

where iV^^^^ iE+, E-) is the distribution of data and N^l^^ is the total number of events in 
that bin. The values of ratios r'*) obtained must be corrected since the MC cannot simulate 
hadronic interactions accurately enough. By introducing mis-ID probabilities, P-KTr^^ti and 
P^l,^l^■K^T, the r value for each bin (the bin number i is omitted) can be written as: 

N + N P — N P 

^ " N Tn ' ^ ' 

where A'^^^ (^/i/i) is the number of true tt+tt^ (//+^^) pairs in that bin. We assume that 
P-rr-K^fifi and Pfi^^TTTt are independent of VF. Applying the /z/tt separation method mentioned 
above to a sample of data events positively identified as muons by the KLM information 
in the higher energy region, we find that P^^-tTTTr is statistically consistent with zero. The 
values of Ptttt^pai in each | cos 9* \ bin are determined such that the ratio of the data and 
MC for fJ,^ ^J,^ pairs, which is ideally one, gives a straight line in the W spectrum. The 
values of Ptttt^/^/j vary between 0.08 to 0.13 in | cos 61*1 bins. Because they are determined 
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for each bin of | cos9*\, the bin-by-bin variation of systematic errors is rather large in the 
angular distribution. After subtracting events, a total of 6.4 x 10® events remains in 

the region of 0.8 GeV/c^ < < 1.5 GeV/c^ and | cos6'*| < 0.6. 

V. DETECTION AND TRIGGER EFFICIENCY 

The detection (trigger) efRciencies, edct (etrg) are estimated from a MC simulation. Events 
of the process 77 tt+tt" are generated using TREPS The detection efficiency is 

calculated from the MC simulation as the ratio of the number of detected and generated 
events in each bin of W (with the bin width, 5 MeV/c^) and |cos6'*| (0.05). The MC 
statistics are high enough and do not contribute to systematic errors. 




FIG. 3: The combined detection and trigger efRciencies as a function of W and | cosS*|. 

The trigger efficiencies are estimated with the same binning using the trigger simulator. 
Since the trigger simulator does not simulate triggers very accurately, particularly in the low 
energy region, the efficiency values have to be corrected. We calculate the correction factors 
by comparing the number of e~^e~ — > e~^e~e~^e~ events in data and MC that are triggered 
by the two-track trigger. The resulting factors steeply rise from 0.5 at W = 0.8 GeV/c^ to 
0.8 at W — 1 GeV/c^ and then increase gradually for higher W. The combined detection 
and trigger efficiencies are shown in Fig. [3l The muon-background subtraction and all the 
correction factors are applied using smooth functions obtained by parameterizing the results 
of bin-by-bin analyses. 
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VI. CROSS SECTIONS 

In this section, we derive differential and total cross sections and evaluate systematic 
errors. 



A. Differential Cross Sections 

Differential cross sections for 77 tt+tt^ are evaluated by using the following relation: 

A\cos0*\ ~ etrg-edofAW-A\cose*\-j§- f Ldt ' 

where AN^+^-^^+^-t^+t^- is the number of events in a W-\ 008 6**1 bin, is the two-photon 
luminosity function and / Ldt = 85.9 fb~^ is the integrated luminosity. Here the W and 
I cos 6** I bin sizes are also chosen to be 5 MeV/c^ and 0.05, respectively. Background from 
r7'(958) p°7 7r+7r~7 is subtracted, a detailed account of which is given in Appendix [Al 
The contribution of the background to the cross section is about 5% at 0.8 GeV/c^ and dies 
away quickly to zero above 0.9 GeV/c^. Other backgrounds are neghgible. 

Differential cross sections da/d\cos6*\ are obtained using Eq. ^ for |cos0*| from to 
0.6 and for W from 0.8 GeV/c^ up to 1.5 GeV/c^. The resulting differential cross sections 
are shown in Fig. |4l In order to present the cross sections more quantitatively, some rep- 
resentative ones are also plotted in Fig. [H Both statistical and point-by-point errors are 
shown. The latter come from the /z/tt separation method and trigger efficiency corrections 
we employ as explained in Sections IIVI and |Vl A point-by-point systematic error is taken 
to be one half of the difference between the corrections in neighboring bins. The result 
shows some apparent systematic structure in the region 0.45 < |cos6'*| < 0.6, particularly 
for W < 1.1 GeV/c^. As shown below (see Eq. ([4])), the differential cross sections in this W 
region can be described by a second order polynomial in | 008 6**1^. Thus, such structures are 
not considered to be real; either the last points are too low or earher points are too high. 
However, we have not identified the cause of the measurement bias. 




FIG. 4: Differential cross sections (dcr/dj cos | (nb)). (The W^-axis is reversed compared to that 
of Fig. Oso as to allow a clearer view of the region in W above the /2(1270) resonance.) 
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FIG. 5: Representative differential cross sections (dcr/d| cos6'*| (nb)). In the figures, 0.80 GeV 
means a bin of 0.800 GeV/c^ < W < 0.805 GeV/c^, etc., and the dashed lines indicate the upper 
and lower overall systematic errors. The two short horizontal bars indicate the statistical errors 
while the vertical ones include point-by-point systematic errors. 



In this W region, J > 2 partial waves (the next one is J = 4) may be neglected so that 
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FIG. 6: The differential cross section of 77 tt+tt" at the /2(1270) mass (1.275 GeV/c^). Thick 
vertical bars show statistical errors and thin ones include point-by-point errors. The dotted lines 
indicate the overall systematic errors. The solid line shows the angular dependence of Y2 and the 
dashed one shows that of Y2 (both normalized at cosO* = 0). 

only S and D waves are to be considered. The differential cross section can be expressed as: 

^^^^^^+^-)^\SY,' + DoY,'\' + \D,Yi\' , (4) 

vifhere Do (^2) denotes the helicity (2) component of the D wave and FJ" are the spherical 
harmonics: 

Since \Y2 \ is not independent of Y^ and Y^ (i.e. jFj^l = (V^Y^ — Y^)/V6), partial waves 
cannot be separated from the differential cross sections alone; additional inputs or assump- 
tions are needed. The general trend of the angular distribution as a function of W is as 
follows. The angular distribution below 1 GeV/c^ is rather flat for | cos 6**1 < 0.4, indicating 
that the S wave fraction is significant (Fig. [5]). In the region above 1 GeV/c^, the angular 
dependence becomes steeper as W increases and is the steepest around the /2(1270) mass. 
Such behavior is typical of D wave dominance. Theoretically, the helicity=2 wave {D2) is 
expected to be dominant [l^. This is supported by Fig. [6l where the angular dependence 
of 1^2^12 and (Y^)^ is plotted at the /2(1270) mass. 

B. Total Cross Section 

The total cross section is then obtained by integrating the differential cross sections over 
I cos6'*| up to 0.6 and is shown in Fig. [7] together with the results of some past experiments. 
A clear peak corresponding to the /o(980) meson is visible, indicating that the peak solution 
of the BP analysis is preferred. Systematic errors for the total cross section are summarized 
in Table HI They are dominated by the uncertainty in the /x/tt separation and that of 
the trigger efficiency. Systematic errors arising from the ^/tt separation are estimated by 
changing the value of P-mv-,^ni in the allowable range in each angular bin. Since events 
are well identified by the KLM for W > 1.& GeV, the allowable range is determined in this 
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FIG. 7: The total cross section of 77 -> tt+tt" between 0.8 and 1.5 GeV/c^ for ] cos 6**1 < 0.6. The 
Belle data are represented by crosses with statistical error bars, the Mark II data are squares, the 
CELLO data are the open triangles and the VENUS data are the filled triangles. Dashed lines 
indicate upper and lower systematic uncertainties for the Belle data. Numerical values are listed in 
Appendix [B] We do not show systematic errors for the other experiments; they are of similar size 
or larger. 



region. These well identified events are also used in estimating systematic errors of the 

trigger efficiency. Comparing data and MC for events in the region W > 1.6 GeV and 

extrapolating linearly downward, the systematic errors are found to be 4% at = 1.5 GeV 
and 10% at 14^ = 0.8 GeV. The total systematic error is obtained by summing the systematic 
errors in quadrature and is also shown in Fig. [71 Our results are in good agreement with 
past experiments except for the /2 (1270) mass peak region, where our data points are about 
10 to 15% larger, but still within the systematic errors. 



TABLE I: Summary of systematic errors for the 77 — > tv^tt cross section. A range is shown when 
the uncertainty has W dependence. 



Parameter 


Syst. error (%) 


Tracking efficiency 


2.4 


Trigger efficiency 


4-10 


K/ TT-separation 


- 1 


/^/TT-separation 


5-7 


Luminosity function 


5 


Integrated luminosity 


1.4 


Total 


11.1 - 12.3 
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VII. FITS TO THE CROSS SECTIONS 



The results of the cross section measurements can be used to obtain the parameters of 
/o(980) and /2(1270) resonances, and to search for other states decaying into tt+tt". Some 
of us plan to perform a full amplitude analysis in the near future using the present data 
including the differential cross sections along with published cross section data of the past. 
Thus, we restrict our analysis to a simple level in this paper. In this section, we summarize 
the measurement of the parameters of the /o(980) discussed in a separate paper fl2|, perform 
a simple fit for the /2(1270) resonance as a consistency check, and search for P- and CP- 
violating decay of the 77'(958) meson into a tt+tt^ pair. 



A. The fo (980) Resonance 



We have to take into account the effect of the KK channel that opens within the fo (980) 
mass region. The fitting function for the scalar resonance /o(980) is parameterized as follows: 



,BG 



BG 



,BG 



(6) 



where is the amplitude of the /o(980) meson Jl^l, which interferes with the helicity- 
0-background amplitude yfTo ^ith a relative phase tp, and cr^^ is the total background 
cross section. The amplitude can be written as 



^/4:8^ 



1 



W Wtt 

where the factor 4.8 includes the fiducial angular acceptance |cos0*| < 0.6, f3x 



(7) 



1 



is the velocity of the particle X with mass mx in the two-body final state XX, 



and gfoxx is related to the partial width of the /o(980) meson via Txxifo) 
The factor Dfg is given as follows jUj: 
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fo 



Df,{W) = ml ~W' + ml" {mf„) - n^" (W) + 3?n^ (m/J - (W) , 
where for X — n or K, (m/„ ) is the real part of H;^ [Tnfa), which is given by: 



n$(w^) = 



XX 



TT 1+Px 



(8) 



(9) 



The phase factor is real in the region W > 2mK and becomes imaginary for W < 2mK- 
The mass difference between K"^ and {K^) is included by taking (3k — ■^{(3k± + Pkq )■ 
The results of the fit (shown in Fig. [Hand in TablelTTl) are discussed in a separate paper [ll|. 
For completeness we report here the parameters of the /o (980) meson obtained from the fit. 



= 985.6 



(stat) Ii.6(syst) MeV/c^ 



(/o) = 34.2 tin (stat) tli (syst) MeV 



(/o) = 205 til (stat) (syst) eV. 



310 tl% (stat) eV, and the value 



The two-photon width given by the PDG Q is F^-y (/o) 
found by BP is 280t^3o Our value of the two-photon width is consistent with them 
within errors. 
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FIG. 8: Results of the fit: (a) the total cross section (the solid curve) (b) contributions of the 
resonance ((7(77 /o(980) ^ 7r~'"7r")) (solid line) and the interference (dashed). The cross section 
of cr(77 ^ /o(980) ^ K'^K-) is also shown (dotted) [3. 



TABLE 11: Fitted parameters of the /o (980) region to Eq.©. 



parameter 


value 


error 

stat 


syst 


m/o (MeV/c^) 
Stttt (GeV) 
r77(/o) (eV) 

(nb) 
f (rad) 


985.6 

1.33 

205 

3.7 

1.74 


+1.2 
-1.5 

+0.27 

-0.23 

+95 

-83 

+ 1.2 

-1.5 

±0.09 


+1.1 
-1.6 

+0.16 

-0.05 

+ 147 

-117 

+4.3 

-3.9 

+0.04 

-0.34 


xV"4f {rid]) 


0.90 (15) 







B. The /a (1270) Region 

Prom the past experiments [3, H, @, 0, S, S it is well known that the position of 
the /2(1270) resonance peak in two-photon production is shifted to lower mass because 
of interference with non- resonant background [22|. In this paper, we give the result of a 
simple fit made as a consistency check in the /2(1270) region. The relativistic Breit-Wigner 
resonance amplitude ^_r(VF) for a spin- J resonance R of mass to_r is given by 



where Fj is the factor coming from the limited solid angle (|cos6'*| < 0.6). Hereafter 
we consider the case J = 2 (the /2(1270) meson). The factor F2 — 0.884 is obtained 
assuming helicity-2 dominance [l^; the angular dependence is assumed to be ■ The 
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energy-dependent total width rtot(W^) is given by 



X 



where X is tt, K, 7, etc. The partial width Txxi^) is parameterized as [23|: 



TxxiW) = TRBiR^XX) 



qx{m%)J D2{qx{m%)rR) ' 



(11) 



(12) 



to 



where Tr is the total width at the resonance mass, qx{W^) = y^W^/A - 
1/(9 + 3a;2 + x'^), and tr is an effective interaction radius that varies from 1 GeV" 
7 GeV~^ in different hadronic reactions [131 ■ For X = tt, K, and 7, the branching fractions 
are OMStoml, 0.046 ±0.004, and (1.41 ±0.13) x IQ-^, respectively 0. For the An and the 
other decay modes, r4^(lF) = TrB[R ~f An)^^ is used instead of Eq. Ill2p . 
The fitting function for the /2(1270) region is taken to be as follows: 



\W)e 



i4>2 



W 



1 GeV/c2 



Co + ClW + C2W'^ 



(13) 



where the contribution other than that of the /2(1270) resonance is subdivided into the 
interfering part (helicity=2) and the non-interfering part (helicity=0). The fit region is 
chosen to be ±rtot around the /2 mass, i.e. 1.090 GeV/c^ < W < 1.461 GeV/c^. The 
parameters of the /2(1270) meson are fixed to the values from the PDG: the branching 
fractions as listed above, = 1275.4 ± 1.1 MeV/c^ and = 185.2t2;5 MeV [2], and the 
parameter is fioated. 

The result of the fit is shown in Fig. [9] and the obtained parameters are summarized 
in Table IIIH where errors shown are statistical only. Since a good fit is obtained with 
C2 = 0, we omit C2. A fit without the non-interfering background gives much worse results 
as summarized in Table Hill We conclude that the consistency check is satisfactory. 
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FIG. 9: Results of the fit of the /2(1270) region to Eq. The parameters of the /2(1270) mesons 

are fixed to the values by the PDG (with the helicity=0 background). 
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TABLE III: Fitted parameters for the /2 (1270) region to Eq. QS]). The parameters of the /2(1270) 
meson are fixed to the values by the PDG. Errors shown are statistical only. 



parameter 


with hel.=0 bgd 


without hel.=0 bgd 


TR (GeV-1) 


3.62 ± 0.03 


2.84 ±0.04 


60 (V^) 


5.54 ±0.02 


7.70 ±0.05 




0.61 ±0.05 


2.00 ±0.02 


h (deg.) 


28.7±0.2 


22.8 ±0.1 


Co (nb) 


18.2 ±3.7 


(fixed) 


ci (nb/GeV/c2) 


-1.7±2.8 


(fixed) 


C2 (nb/(GeV/c2)2) 


(fixed) 


(fixed) 


xV"4f {ndf) 


1.1 (68) 


2.1 (70) 



C. The tt'^tx branching fraction of the rf' (958) meson 

The rj' (958) meson is a pseudoscalar meson and, thus, its coupUng to tttt violates P and 
CP. The present upper Umit for the tt+tt^ branching fraction B{ri' — > tt+tt") is 2% (25| . 
The high statistics data of Belle allow for a more sensitive search. The 77' (958) meson has a 
small width of T^' = 0.202 ± 0.016 MeV and a mass of m^- = 957.78 ± 0.14 MeV/c^. Thus 
its contribution to the W spectrum can be represented by a Gaussian function: 

fAW)dW = -^exp f- ^^-T"'^' ) dW , (14) 

where S^i is the parameter to be determined, and aw = 2.0 ± 0.2 MeV/c^ is the mass 
resolution determined from MC. The total cross section in the region 0.92 GeV/c^ < < 
0.98 GeV/c^ is fitted with a second-order polynomial plus Eq. (fT4l) . The finite bin-size 
effect is taken into account by integrating the Gaussian over each bin. The result is S",,' = 
—27 ± 16 nb • MeV. The systematic error is found to be negligible, which is estimated by 
constraining the 77' mass and the mass resolution within one standard deviation and by 
changing the fitting region. 

The parameter S^i can be related to the tt+tt^ branching fraction S(r/ —* tt+tt^) as 
follows. The cross section formula to be used is the same as Eq. ([6]) except for replacing the 
amplitude T^'' by Eq. JTO]) with J = 0: 

4.87rm^> ^V.^,Y ^^B[r]' ^ 7r+^ ^ \.%^V ^^B{j]' tt+^-J 

where T^^ — 4.30 ± 0.15 keV is the two-photon width of the t]' meson, and the latter 
equation is obtained in a narrow width approximation. Taking into account an interference 
effect (for the most conservative case) and using the relation dW / {{W - m,,f + Tl, / A) ~ 
27r/r,j/, we obtain: 



^^^^ ^ L2^ / 2.T,,B{i^.+.^) ^ ,i,^y2vraBGr,,r,,S(v^.+.-)) , (16) 

where ctq'^ is the cross section of the continuum 77 tt+tt^ component whose amplitude 
interferes with the P- and CP- violating 77' decay, and ip' is the phase angle and sin 1^9' = — 1 
gives the most conservative upper limit of S,' . 

To obtain the upper Hmit for S,' at 90% confidence level (C.L.), we have to consider 
two physically possible cases: S,' is negative or positive, depending on the presence or 
absence of an interference effect between amplitudes of opposite P and CP. As the reaction 
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FIG. 10: Fit result of a Gaussian + second-order polynomial in the r;'(958) region. The solid 
(dotted) line is the 90% C.L. upper limit without (with) interference. 

77 ~^ v' would take place via a P wave, while only even orbital angular momentum waves 
can contribute to the ordinary 77 tt+tt" process, it is unlikely that these two processes 
would interfere. In that case, ctq'^ — and Sjj' is non negative. On the other hand, if 
interference is present, the lowest boundary of S,f is — Tro-Q^r^'/^ — —30 nb-MeV, where 
ctq ^ = 93.5 nb is used, i.e. the largest possible value of ctq'^ that gives the most conservative 
limit. 

We first obtain S^? , the 90% C.L. upper limit of S*,,' from the following relation: 

£-;exp(-^!(M).s,, .0^9£»p(-^!<M).5,, , ,17) 

77' r]' 

where X^i^i]') is the from the fit with a fixed S^f and S*™™ is the lower physical boundary 
of Sjf'. In the presence (absence) of the interference effect, the limit is determined to be 
5"^]? < -2.0 (S"^? < 14.4) nb-MeV. The results are shown in Fig. [H We obtain the upper 
limit of B{r]' — > tt+tt^) taking the errors of F^^ and B^.y (= F-y^/F^') into account. Namely, 
we calculate the contribution to the uncertainty in Srj' that arises from these parameters 
and combine it with the statistical error of Srj', 16 nb • MeV, reevaluate S"^?, and then 
translate it into a limit for B{r]' — > tt+tt^). In the case of no interference, we obtain 
B{ri' tt+tt^) < 3.3 x lO"''. In the other extreme case of maximum interference, we use 
cr^° = 93.5 nb, and the limit is B{r]' tt+tt") < 2.8 x lO'^ at 90% C.L. The errors in F^-^ 
and Bj^ are also included but they lead to a negligible change in the upper Hmits. 

Vm. SUMMARY AND CONCLUSION 

In summary, we have performed a high statistics measurement of the 77 tt+tt^ cross 
sections in the tt+tt" invariant mass region 0.80 GeV/c^ <W<1.5 GeV/c^ for | cos 6**1 < 0.6 
with the Belle detector at the KEKB e+e" collider. The total cross section is measured in 
fine bins of W {AW = 5 MeV/c^) and differential cross sections are given in bins of AW = 
5 MeV/c^ and A\ cos 6**1 = 0.05. We have observed a significant peak corresponding to the 
/o (980) resonance. Our data clearly select the peak solution of the Boglione-Pennington 
amplitude analysis [ll|. The total cross section is fitted to obtain the parameters of the 
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/o(980) meson [14] and to check consistency in the /2(1270) region. For a P- and CP- 
violating decay of the 77' (958) meson, we set an upper limit without (with) interference 
between opposite P and CP amplitudes B{tj'{958) n+n") < 3.3 x 10"'' (< 2.9 x 10"^) 
at 90% C.L. thereby significantly improving the previous Hmit of 0.02 The angular 
dependence of the differential cross sections is consistent with the presence of a significant 
S wave fraction for W < 1 GeV/c^ and with the dominance of the D wave in the /2(1270) 
mass region. 
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APPENDIX A: BACKGROUND SUBTRACTION 

In this section, we describe in detail how the background from 77' 7r+7r~7 is subtracted 
in the cross section determinations described in Section IVII and Appendix [BJ Note that 
beam-gas background that was important in past experiments is completely negligible at 
B-factories because of the very high luminosity. The dominant physics background is due 
to 77'(958) production and its subsequent decay into /5°7; the photon energy at the nominal 
/o" mass is 0.14 GeV, and a significant fraction of tt+tt" pairs satisfy selection requirements 
such as I |< 0.1 GeV/c. The other physics backgrounds are neghgibly small. 

The background from 77' 7r+7r~7 can be estimated with MC at the four- vector level 
without doing a full detector simulation. This is because the subtractions are applied to 
cross sections where efficiency corrections of the detector and trigger are already included 
and because the bin-sizes used are large enough that the effect of finite detector resolution 
is neghgible. The process e~^e~ e~^e~r]' with 77' p^j 7r"'"7r~7 is simulated using 
TREPS [l3|- The final state n^n^j is generated with a matrix element incorporating the 
dipole transition feature of the 7/ p^j decay and the p pole (25l . [2^ : 

piElm? , _ sin^ 6 

where is the pion momentum, is the photon energy and 6 is the angle between one of 
the pions and the photon, all evaluated in the di-pion rest system, and the denominator is 
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the p pole. The mass dependence of the p width is parameterized as [22 



r(m) 



ro , 



(A2) 



where po is the p^r at to = TOp and Fq is the nominal p width. Generated tt+tt" pairs are 
subjected to the cut | X]Pt l< ^-^ GeV/c and then accumulated into bins oiW - |cos0*| 
with the same bin size as that of the differential cross sections. The obtained distribution 
is related to the cross section of 77 ^ X by using [18]: 



7 r* 

da{e+e- ^ e+e'X) = dW—iW)a{-f-f X) 



(A3) 



The resulting background cross sections to be subtracted are shown in Figs. 11(a) and (b). 
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FIG. 11: Background total and differential cross sections (77 -rj' ^ p^'j tv^tt 7) 
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APPENDIX B: TOTAL CROSS SECTION 

In this appendix, we list the measured values of the total cross section for the process 
77 tt+tt" integrated over the angular region | cos9*\ < 0.6 in the range of 0.8 GeV/c^ < 
W < 1.5 GeV/c^ in steps of AW = 0.005 GeV/c^. At each energy, the third number in the 
table is the statistical and the fourth is the overall systematic error. 
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